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1. Introduction 

This talk summarizes some recent theoretical developments in high energy QCD, 
in particular those that deal with very high energy heavy ion collisions. Compared 
to the early days of studies of heavy ion collisions, the developments have been 
significant though outstanding problems remain. One such development is a QCD 
effective field theory (EFT) of high energy wavefunctions, the Color Glass Conden- 
sate (CGC) PQ, which provides a unified systematic non-perturbative (albeit weak 
coupling) approach to compute key features of e+A, p+A and A+A collisions. This 
framework has been applied successfully to describe inclusive and diffractive data 
at HERA, inclusive data from fixed target e+A experiments, singly and doubly 
inclusive deuteron-gold data from RHIC, the energy dependence and multiplic- 
ity distributions of p+p collisions up to the highest LHC energies, the centrality 
and multiplicity distributions in A+A collisions, and their energy dependence from 
RHIC to LHC. Long range rapidity correlations in p+p and A+A collisions arise 
naturally in the CGC; a noteworthy recent example is the "ridge" in p+p colli- 
sions which was predicted in this formalism. While these comparisons to data are 
not fully quantitative, the economy and consistency of parameters, and the unified 
description of multi-particle production in deeply inelastic scattering and hadronic 
collisions, make the CGC an attractive candidate for a quantitative description of 
the formation and evolution of quark-gluon matter in heavy ion collisions. 

We shall first briefly discuss the many-body structure of high energy wavefunc- 
tions in the CGC. We next describe how these wavefunctions collide and decohere 
to produce quark-gluon matter called the Glasma. The Glasma, initially described 
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by high occupation number fields, evolves in a highly non-trivial manner to likely 
form a thermalized Quark-Gluon Plasma (QGP). In the final section, we will outline 
a scenario whereby this process occurs, and shall discuss possible implications of 
these results. 



2. Before the collision: hadron/nuclear structure at high energies 

QCD in high energy Regge asymptotics is a dense many-body system of "wee" 
gluons and sea quarks. In the infinite momentum frame, gluons with transverse 
momenta k± < Qs saturate phase space maximally, where Qs(%) is a dynamical 
saturation scale [5] that grows with decreasing fractions x of the longitudinal mo- 
mentum of the hadron carried by the gluons. The properties of saturated gluons 
are described by the CGC EFT: the degrees of freedom arc static color sources 
at large x, coupled to dynamical wee gluon fields at small x. Because of the large 
occupancy of wee gluons, the effective ground state wavefunction is a classical non- 
Abelian Weizsacker- Williams field [3], whose properties at small x, on account of 
Qs Aqcd, can be described in weak coupling. Renormalization group (RG) equa- 
tions, derived from requiring that observables be independent of the separation in 
x between sources and fields, lead to an infinite hierarchy of evolution equations 
in x, for n-point Wilson line correlators averaged over dense color fields in the nu- 
cleus. For a physical observable defined by an average over all static color source 
configurations, 

(0) Y = J[Dp}W Y [p]0[p], (1) 

the gauge invariant weight functional of color sources W Y [p] at rapidity Y — 
ln(j;/a;o) = Id.(xq /x~), where xq is the initial scale for small x evolution, satis- 
fies the JIMWLK equation g] dW Y [p]/dY = UW Y [p]. The energy evolution of the 
observable can then be expressed as 

= (n o) Y . (2) 

The structure of the JIMWLK Hamiltonian % is such that (HO) Y is an object 
distinct from (0) Y ; one obtains an infinite hierarchy of evolution equations for 
operator expectation values (0} Y . Given appropriate initial conditions at large x, 
solutions of this Balitsky- JIMWLK hierarchy [5j 0] allow one to compute a wide 
range of multi-particle final states in deeply inelastic scattering (DIS) and hadronic 
collisions. 

The simplest example is that of the inclusive DIS structure functions F2 and 
Fl , which are proportional to the forward scattering amplitude of a qq "dipole" on 
a nucleus. The forward dipole amplitude (dipole cross section) is given by 



a dip .(x,r T ) = 2 J d 2 b T x <^l--lTrv(b T + ^) Vt ^ 



T, 



(3) 
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where Tt = xr — yr is the transverse size of the dipole, hr — (xt + yr)/2 is the 
impact parameter relative to the hadron and V(xt) = Pexp(ig J dz~ ^ X ya •* ). 
From eq. ([2]), one can show that the expectation value D = (D), with D(xt — yr) = 
jj- Tr V(xt) V (Yt), obeys the Balitsky-JIMWLK equation that relates its energy 
dependence to the expectation value of a four-point operator, 

-^D(x T - y T ) = jf /C xyz (£>(x T - zt) ^(zt - yr) - D(x T - y T )) ,(4) 

where /C xyz = (xr — yr) 2 /( x t — zt) 2 (zt — yr) 2 - In the large A c approximation, 
the expectation value of D 2 factorizes-the resulting closed form expression 

d N c a s f 

— D(x T - yr) = ^ 8 y ^Cxyz [£>(*t - zt) £>(zt - yr) - f(xr - yr)] • (5) 

is known as the Balitsky-Kovchegov (BK) equation [5J [B]. In addition to eq. ([3]), 
the dipole correlator appears in a number of final states in both DIS and hadronic 
scattering; the BK equation for its energy evolution is successfully employed in 
phenomenological applications. 

For less inclusive observables, one encounters the expectation value Q = (Q) of 
the "quadrupole" operator 

Q(y T ,x T ,u T ,v T ) - ^Tr (vHy T )V (x T )V\u T )V (v T )) ■ (6) 

Unlike (DD) in eq. ((4j) , it is not reducible to the product of dipoles even in the large 
N c and large A approximations and is a novel universal correlator in high energy 
QCD [71 IB], interesting both from theoretical and phenomenological perspectives. 
Examples of where this quantity appears include small- a; di-jet production in e+A 
DIS^, quark-antiquark pair production in hadronic collisions [5] and near-side long- 
range rapidity correlations |10) . Another interesting quantity 8 is a six-point corre- 
lator S*6 oc (QD) that is probed in forward di-hadron production d + A — > hi h% X 
in deuteron-gold collisions at RHIC. The dominant underlying QCD process is the 
scattering of a large x\ valence quark from the deuteron off small xi partons in the 
nuclear target, with the emission of a gluon from the valence quark either before 
or after the collision. The RHIC experiments show that the away-side peak in the 
di-hadron correlations is significantly broadened for central collisions at forward 
rapidities [TTJ [12] as predicted in the CGC framework [13] and confirmed in more 
detailed analyses [14l [15] . These analyses however relied on simplified factorization 
assumptions that are not justified. The JIMWLK RG equations for Q [16] and 
Sq |17) have been derived and computed numerically for particular configurations 
of these operators |18) . A remarkable result of these numerical simulations is that 
the initial condition and the RG evolution of these higher point correlators are well 
reproduced by expressions that are functions only of the dipole expectation value 
-D-see ref. [19] for a recent analytical interpretation of these results. 

As this discussion suggests, significant progress has been made in understanding 
the RG evolution of novel multi-point correlations that comprise the structure of 
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hadron wavefunctions at high energies and observables sensitive to these correlators 
in p+A and e+A collisions have been identified. An outstanding problem that is 
unresolved is the impact parameter dependence of high energy evolution. Exclusive 
measurements at an e+A collider provide the best opportunity to shed empirical 
light on this topic 20 . 

3. Creating bulk matter: The Glasma 

Before the collision, the effective ground states of the nuclei are coherent states 
described as classical fields with occupation numbers 0(1/ g 2 ). These fields become 
time dependent immediately after the collision. Inclusive quantities, computed in 
the presence of these time dependent color fields, can be expressed in terms of 
retarded propagators, thereby allowing real time computations with initial data 
at negative infinity |21[ 122] . Other examples of similar non-perturbative dynamics 
include Schwinger's mechanism for electron-positron pair production in strong QED 
fields and Hawking radiation in the vicinity of the event horizon of a black hole. 

The dynamics of the Glasma [29] created in heavy ion collsions is analogous to 
important aspects of early universe cosmology, one being an"Ekpyrotic" picture of 
the production of multiverses through the collision of brancs; another powerful anal- 
ogy is the dynamics of pre-heating that leads to thermalization and hydrodynamics 
of classical fields [24]. In QCD, a classical description of Ekpyrosis is obtained by 
solving Yang-Mills equations with static light front color sources (25] [26] [27] . At this 
classical level, energy dependence is introduced "by hand" via the saturation scale 
Qs- The energy dependence arises at next-to-leading (NLO) order from quantum 
fluctuations (of relative strength g) about the classical background of the nuclei. 
These naively sub-leading contributions are enhanced by logarithms, which give 
asln(xi i 2) ~ 1 for small and have to be resummed to all orders in perturba- 
tion theory. In the strong field regime, there is an additional resummation (gpi^) 71 
at each order in the as expansion. These respectively radiative and multiple scat- 
tering contributions are generated by the JIMWLK RG equation for the weight 
functionals Wy 1{2) [pi(2)] f° r each of the incoming nuclei. First computations of the 
inclusive multiplicities in heavy ion collisions by solving Yang-Mills equations with 
JIMWLK generated initial conditions have been performed recently [28] . 

An important feature of the Glasma is that it generates long-range correlations 
that are localized on scales ~ 1/Qs in the transverse plane [30]. The long range 
correlations evolve with energy and the the rapidity separation between correlated 
gluons [31]. This formalism predicted [32] a near-side "ridge" in high multiplicity 
collisions, an effect that was observed by the CMS collaboration [33]. The system- 
atics of the data are in qualitative agreement with the Glasma predictions [34] . In 
A+A collisions, the observed ridge is a consequence of the boosting of long-range 
correlations in the final state by radial flow [35]. The radial flow provides the near 
side angular collimation. The combination of flux tube structures in the initial state 
and radial flow give a good description of RHIC and LHC data [36] [37] . While it 
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would be interesting to discuss at length these long range correlations, the rest of 
this talk will focus on the space-time evolution of the Glasma. 

Albeit causality requires that nuclei don't communicate before the collision, it 
is by no means assured in a weak coupling treatment; we are able to show formally 
that factorization of the contributions of the weight functionals W[pi,2] to inclusive 
quantities is obtained at leading logs in x accuracy [351 [Ml ED]- These quantum 
modes are boost-invariant p v — modes, where p v is Fourier conjugate to the 
space-time rapidity. In the process of Ekpyrosis, nuclear coherence is lost and p v ^ 
modes are generated. These modes are generically unstable [HI [42l |43j [44] , and grow 
in an expanding system as (as exp(2 v / ^))" with \i ~ Q$, where n denotes the 
order in perturbation theory beyond the classical leading order contribution. These 
"leading instabilities" are comparable to the background field at r ~ 1/Qs and 
have to be resummcd to all orders, leading to qualitatively different behavior. 

After factorization and resummation of leading logs in 2:1.2 and leading instabil- 
ities as exp(2y//xr), the energy-momentum tensor is expressed as 

(T^) LLx+L i nst . = J [D Pl D P2 ] W X M W X2 [p 2 ] J[Va] F [a] T^[A[p 1 ,p 2 ] + a](x) . 

(7) 

The argument A = (A, E) denotes collectively the components of the classical 
fields and their canonically conjugate momenta on the initial proper time surface; 
analytical expressions for these are available at t = + 25, 45 . The initial spectrum 
of fluctuations Fo [a] , Gaussian in the quantum fluctuations a, has a variance given 
by the small fluctuation propagator in the Glasma background field as r — > + . In 
practice, the path integral in a is determined by solving the classical Yang-Mills 
equations repeatedly with the initial conditions at r = + given by 



Ainu. = Ainit. + J dfi K ]c K a£(x) + c* <"(x)J . (8) 

Here A collectively denotes the quantum fields and their canonical conjugate mo- 
menta. The coefficients c K , with K collectively denoting the quantum numbers 
labeling the basis of solutions, are random Gaussian-distributed complex numbers 
given by 

(c k c* k ,) = ^5 kk , , ( Vk ,) = (<c;)=0, (9) 

Explicit expressions for the small fluctuations and their conjugate momenta, denoted 
here by cl^(x) were obtained in ref. |46) . The inner product of these solutions satisfies 
the orthogonality condition (a K | a , ) = N K 8 , with the measure d\x K (a mix of 
integrals and discrete sums) that ensures J dfi K N K 5 KK , = I. 

In ref. [46 , a numerical algorithm was outlined to compute eq. (J7J|, thereby 
describing both Ekpyrosis and inflationary dynamics including essential leading 
quantum fluctuations. (We emphasize the formalism holds for any inclusive quantity 
in heavy-ion collisions including parton energy loss and sphaleron transitions at 
early times.) We can now study how lumpy initial conditions for color charges p\ 2 
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in nuclear wavefunctions at the energy/rapidity of interest transform into the flow 
of matter. The early universe analogy is becoming more robust experimentally with 
a "WMAP-like" spectrum of spatial anisotropies resolved in heavy ion data [27] . 

We carried out an extensive study of the formalism outlined here for a scalar </> 4 
model which, among several QCD-like features, has a spectrum of unstable quantum 
modes which are amplified by resonant interactions with the background field [48) . 
In (f) A (and other scale invariant theories), the amplitude of the field is inversely 
proportional to its period of oscillation. Slightly different amplitudes, correspond- 
ing to different quantum seeds, lead to differing oscillation periods; a stochastic 
average over these leads to decoherence in the evolution. A striking consequence 
is hydrodynamic flow with an ideal equation of state. Eq. ([8]) is a realization of 
Berry's conjecture [49] which is believed to be necessary for thermalization [50] of 
a quantum system. Thermalization and onset of quasi-particle dynamics have been 
studied in the scalar theory for a fixed box [51]; similar studies are feasible in the 
expanding case and eventually in the QCD framework of ref. [46]. A promising de- 
velopment [52] is an attempt to connect the CGC power counting approach here 
with a well developed 2PI formalism [53j [54] which has to potential to allow one 
to follow the evolution of the system in heavy ion collisions until thermalization. 
In the following section, we will describe a complementary approach based on the 
Boltzmann equation that describes how the highly occupied initial state approaches 
thermalization. 

4. Space-time evolution of strongly interacting matter: from 
Glasma to Quark-Gluon Plasma 

The classical Yang-Mills energy-momentum tensor in the Glasma is of the form 
T^ v = diag (e, e, e, — e), and therefore has a negative longitudinal pressure at very 
early times [55] . The unstable quantum fluctuations we discussed above can however 
rapidly isotropize the energy-momentum tensor on very short time scales of order 
1/Qs- An outstanding issue then is whether the system can resist the tendency of 
to fall out of isotropy again due to its rapid expansion into the vacuum. Another 
outstanding issue is whether the system can equilibrate to generate a Bose-Einstein 
distribution on the short time scales available. The "bottom-up" scenario [56] out- 
lined a systematic way whereby hard elastic and inelastic collisions of produced 
gluons would drive the system towards equilibration. In ref. [57], it was demon- 
strated that anisotropy driven instabilities can significantly alter the bottom-up 
picture. Here, we will discuss an alternative scenario |58j where early time instabil- 
ities, as discussed previously, have already played a role, generating a gluon density 
in the Glasma that is parametrically large compared to the value it should have 
in a system in thermal equilibrium with the same energy density. In such systems, 
the excess of gluons can be diluted by dynamical generation of a Bose-Einstein con- 
densate, corresponding to a large occupation of the zero momentum mode, and/or 
by inelastic processes that in the long run tend to tame the particle excess. Until 
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the latter begin to dominate, a transient Bose-Einstein condensate can exist, as in 
superfluids, with interesting consequences. 

We will assume that after times ~ 1/Q S , the Glasma can be described by color 
singlet distributions for both the particle content and the condensate. We will study 
in kinetic approach the role of collisions in driving the overpopulated system to 
equilibration. We will assume that, at all times t > 1/Qs, the distribution function 
takes the form 

f(p) ~ — forp< A s , f(p) ~ - — - for A s < p < A, f{p) ~ for A < p. 
a s a s w p 

(10) 

At t ~ 1/Qs; both scales A s and A coincide with Q s . As time progresses, the 
two scales separate, with A s decreasing quickly, and A evolving much more slowly. 
Thermalization is reached when A s /A ~ ct s , at which point, /(A) becomes of order 
unity. 

A more precise definition of these scales is obtained from the collision integral of 
the Boltzmann equation in the small angle approximation, assuming 2 — > 2 elastic 
scattering and isotropy, 



dj_ 
dt 



coll P 2 



dp{p 



,2 



I") 



The fixed point solution of this equation is a Bose-Einstein distribution with tem- 
perature T = A s /a s . The scales A s and A are obtained from the integrals 

AA S _ f°° , 2 df A A 



/ d PP±> —T = / d PP 2 f(l + f). (12) 
a s J dp a{ J 

Remarkably, in the regime where / S> 1 (/ ~ l/a s )i all dependence on a s drops 
from the collision integral. Taking moments of the collision integral with arbitrary 
powers of p, the typical collision time is given by i sca t = -£s ■ The scattering time 
is itself a function of time and simple analysis of moments of the kinetic equation 
suggests that i sca t ~ t. For a fixed box, the energy conservation conditioro gives 
A S A 3 ~ constant. With these two conditions, one obtains the temporal evolution of 
the two scales to be 

3 1 

A 8 ~Q S (^)%A~Q 5 (£)\ (13) 

The number density of gluons n g ~ A 2 A S decreases as ~ (to/t) 1 / 7 . The Debye mass 
m = ^/A S A decreases slowly in time as ~ Qsito/t) 1 ^"' ' ■ The thermalization time, 

determined from A s ~ a s A, is t t h ~ ^- (^^j 4 • Because the scale A increases with 

time, so does the entropy density s ~ A 3 ~ Q 3 (t/to) 3 ^ 7 . When t — t t h, s ~ Q 3 /a^ 4 , 



a The gluon number is not conserved because of the B-E condensate and/or inelastic number 
changing processes. We are also assuming here that the energy density of the condensate is small, 
as argued in ref. 58_ . 
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which is the equilibrium entropy ~ T 3 . Finally, becase quarks have a phase space 
density of order 1 up to the scale A, the number density of quarks is n quar k s ~ A 3 . 
At t ~ 1/Qs) this is suppressed compared to the number of gluons n g ~ A s A 2 /a s 
by l/a s . They become of the same order when A s ~ a s A, the thcrmalization time, 
and cannot be ignored at this time. 

We now consider the effect of longitudinal expansion, which can be included by 
adding a drift term to the left hand side of the kinetic equation. By integrating over 
momentum the kinetic equation multiplied by the energy one obtains 

3 t e+^±A = 0, (14) 

where e is the energy density, e = J ^J^ 3 w p f p , and P L the longitudinal pressure, 

Pl = J Sfy f-/ P ' The effect of longitudinal expansion is observed by parame- 
terizing the longitudinal pressure in terms of energy density as P L = S e , where 
the multiplicative factor 8 can be in the range [0,1/3]; 6 — corresponds to the 
free streaming limit and S = 1/3 to that of to ideal hydrodynamic expansion after 
isotropization. The assumption that 6 is independent of time is a strong one which 
we make to understand how collisions redistribute momenta and thereby generate 
the shape of a thermal distribution for an expanding system. 

Within this framework, e g (t) ~ e(t ) (^) + . Our previous estimate of the col- 
lision time is unchanged and can be combined with this expression to yield, 

/+ \ (4+<5)/7 \ (l+2«5)/7 

A s ~Qs(JJ , A~Q S (JJ . (15) 

From these, one easily obtains the estimates of the gluon density, and of the Debye 
mass, 

n 3 / f \ (6+5<5)/7 / x (5+M)/7 



m^Qi[j) (16) 
The thermalization time, obtained as before from the condition A s = a B A, is given 

7 

by ~ (ri - ) 3 S • Comparing this result to the fixed box result, we see that 

the expansion has the effect of delaying thermalization. (Formally, one recovers the 
fixed box result by setting S = — 1, which corresponds to constant energy density). 

For moderate values of the anisotropy (more precisely for 6 > 1/5), a condensate 
can form, with number and energy density given by 

(5-ll(5)/14 



Qt (to 
t 



^'7 



N (-l+55)/7' 



(17) 



The energy density of the condensate is subleading and decreases with increasing t. 
We should also note, as shown in ref. [5S] , that in a relaxation time approximation 
of the collision integral, solutions exist that correspond to a fixed time independent 
anisotropy. We shall now comment on the role of inelastic 2 -H> 3 processes [55] , 
One can show that the inelastic scattering rate is parametrically the same order 



9 



iineiastic ~ j2 ~ ^scat- Because the 2f) 2 transport equation for the overoccupied 
distribution (oc l/w p ) provides a powerful source term for the Bose- Einstein con- 
densate, it can exist as a transient state perhaps until close to thermalization when 
presumably inelastic scatterings become dominant. Finally, we should note that the 
result of ref. [55] for a fixed box coincides with that of ref. [50] for (among many 
scenarios considered) what we consider to be the physical scenario corresponding to 
the Glasma. For an expanding Glasma, the emphasis in ref. |61) is on plasma insta- 
bilities, while in our case strong elastic scattering is the dominant mechanism. Our 
results for thermalization do not require Bose-Einstein condensation; the possible 
creation of a transient condensate is plausible and may have interesting phcnomcno- 
logical consequences. We anticipate that future non-perturbative non-equilibrium 
simulations outlined in ref. [46] (and possible refinements thereoff [52]) will help 
clarify this issue. 
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